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We have performed a one-loop calculation of the width of the rare decay H — > ZZ^j in the 
standard model for Higgs boson masses 190 GeV < mn < 250 GeV. We find that the most dominant 
helicity combinations for the Z bosons and the photon is when one of the Z bosons is longitudinally 
polarized and the other Z boson and the photon have the same helicity. A comparison of the decay 
width T(H —* ZZj) to those of H — > 77 and H —* 7Z shows that the ratios of the decay widths 
are F(H -> ZZ-y)/T(H -» 77) ~ T(H -> ZZj)/F{H -> jZ) < 10~ 7 . 



PACS numbers: 13.15.+g, 14.60.Lm, 14.70.Bh, 95.30.Cq 



1. INTRODUCTION 

In the standard model, the lowest order contribution 
to the decay H — > ZZj takes place at the one-loop level. 
The Feynman diagrams for this process are similar to 
those of the rare decay process H — > 77^0- The decays 
H — > 77 and H — ► 7Z which also occur at the one-loop 
level, and H — * ZZ which occurs at the tree level, dom- 
inate the decay process H — > ZZ^f by several orders of 
magnitude 0. However, by imposing kinematic cuts on 
the Z bosons and the photon in the decay H — > ZZ"/, we 
may exclude contributions of the back-to-back Z bosons, 
and distinguish the Z bosons of the H — > ZZ7 decay 
from those of the H — > decay. Although very rare, 
the decay i? — > ZZ7 is sensitive to the top quark cou- 
plings, and therefore, a signal in this channel is an evi- 
dence for a deviation from the standard model top quark 
couplings. 

In the next section, we outine our calculations of the 
decay width, the Z boson invariant mass decay distribu- 
tion, and the photon's energy decay distribution. This is 
followed by a summary and conclusions. 



—ig^u/ik 2 —m 2 ), where k and m are the momentum and 
the mass of a gauge boson, respectively. 

The resulting Feynman diagrams are similar to those 
of the rare decay process H — » 77Z and its crossed 
channel scattering process 77 — > ZH 4J. Some repre- 
sentative diagrams are drawn |5( in the Fig.^ In Ref. 





FIG. 1: Some representative Feynman diagrams for the decay 
process H — > ZZ7 are shown. For the charged fermion / in the 
Figs.^a) andnTb - ). we include only the top quark. In the Fig.HTc'). 
all charged fermions / of the third generation are included. 



DECAY WIDTH CALCULATION 



In order to decrease the number of the Feynman dia- 
grams needed to calculate the amplitudes for the process 
H — > ZZ~f, we use the generalized non-linear gauge fix- 
ing condition introduced in the Ref. |3( ■ The parameters 
contained in the gauge fixing terms do not affect observ- 
ables and may be chosen for convenience. We choose the 
values of these parameters in such a way as to eliminate 
certain vertices and thereby minimize the number of dia- 
grams to be calculated. To further simplify the structure 
of the resulting tensor integrals and to increase their nu- 
merical stability, we use the 't Hooft-Feynman gauge in 
which the propagators for the gauge bosons have the form 



[j] it was noted that the charge conjugation properties 
of the gauge boson couplings exclude any contribution 
from W bosons in the loops and only charged fermions 
in the loops contribute. In the Fig.^a), the coupling 
of the Higgs boson to the fermion is proportional to the 
fermion mass and we include only the top quark contri- 
bution. The same is true in the Fig.^b) for the coupling 
of the neutral Goldstone boson G° and the fermion, and 
again we include only the top quark contribution. 

The contribution of the diagram of the Fig^c), apart 
from the dependence on the mass of the fermion in the 
loop, is proportional to N[ Qfg^gy- Here, N[ is the 
number of fermion colors {N[ = 1 for a lepton, N[ — 3 
for a quark) , Qf is the fermion electric charge in units of 
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the proton charge, g* A = r/ and g{? — T/ — 2Q f sin 2 Ow 
are the axial- vector and the vector coupling constants, 
respectively. Tg is the third component of the weak 
isospin, and 9\y is the weak mixing angle. As far as 
its dependence on the fermion mass is concerned, the 
contribution of the diagram of the Fig^c) consists of 
two parts. One depends on the mass of the charged 
fermion in the loop (it vanishes for a massless fermion), 
and the other is independent of the fermion mass. The 
latter gives an anomalous contribution [||. However, it 
is clear that the inclusion of all charged fermions of a 
given generation will cancel this anomalous contribution 

since Y^f N% Qf g A 9v = 0- Furthermore, if all of the 
members of a particular generation had the same mass, 
the total contribution of that generation would vanish. 
We included only the charged fermions of the third gen- 
eration in the evaluation of the diagram of the Fig.^c) . 
(The other two generations give negligible contributions.) 

The process H — > ZZ7 will be dominated by the 
decays H — ► 77, H — > 7Z, and especially H — > ZZ 
that occurs at tree level. In order to facilitate the dis- 
crimination of H — > ZZ^f from these dominant decay 
modes and account for some of the possible experimental 
limitations, we imposed cuts on the following Z boson 
and photon kinematic variables: \pz\, \pz'\, \Py\> m zz'i 
m %z> m yz>' °zz', 6yz, and jZ '- Here, p z , p z >, and 
p 1 are the 3-momenta of the Z bosons and photon, re- 
spectively, in the center of mass of the Higgs boson, and 
Qzz'i 6-yZ, and Q^z' are the various angles between the 
3-momenta, pz, pz', and p 7 . The invariant mass vari- 
ables are m 2 zz , = (p z + Pz>) 2 , m^ z = (p 1 + p z ) 2 , and 
m^ z , = (p 7 + pz 1 ) 2 , where pz, pz>, and p 1 are the 4- 
momenta of the Z bosons and the photon. 

We choose the following set of cuts for the numerical 
calculations of the decay width: 

|Pz|cut = |PZ'|cut = l^ylcut = |p|cut, (1) 

(razz')cut = (myz)cut = (mjz')cut = "i cut , (2) 

{OzZ')cut = (#7z)cut = = #cut ■ (3) 

These kinematic cuts facilitate the experimental tagging 
of the Z bosons and photon. They provide minimum 
opening angles between the Z bosons and the photon, 
exclude contributions of the back-to-back Z bosons and 
photons, exclude soft photons, and also improve the nu- 
merical stability of the calculations. The cuts help dis- 
criminate the non-back-to-back Z bosons and the photon 
of the decay H — > ZZj from the back-to-back Z bosons 
and photons in the decays H — » ZZ, H — » jZ, and 
H — ► 77. In principle, all the Z bosons and photons of 
the decays H -> ZZ7, H -> ZZ, H -> 7Z, and H -> 77 
can be identified. However, resolving ambiguities in the 
reconstruction of the Z bosons can be problematic, es- 
pecially when decays into neutrino-antineutrino pairs are 
involved. 

For the calculation of the decay width, and its distri- 
butions with respect to the invariant mass and energy, we 



have used a semiautomatic |8( method. In this method, 
we used the FeynArts package @ to generate the helic- 
ity amplitudes in terms of the tensor loop i nteg rals and 
then, using the Passarino-Veltman method [l(j, we ex- 
pressed these tensor integrals in terms of scalar integrals 
and evaluated them numerically |ll| . 

The result of the calculation for the decay width 
T(H — > ZZ^f) as function of the Higgs boson mass rag 
is shown in the Fig.[2 For comparison, in this figure we 
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FIG. 2: (Color online) The decay widths as function of mn for 
several decay modes of the Higgs boson are shown. The solid line 
is T(H — > ZZj), the dashed line is F(H — > 77), the dotted line 
is r(_ff -» ■yZ), and the dotdashed line is V(H -* ZZ). The cuts 
imposed on F(H — > ZZy) are |p| C ut = 5 GeV, m cu t = 10 GeV, and 
6 cut = tt/24. 

also included the decay widths T(H — > ZZ), T(H — > 77), 
and T(H — > jZ) It is clear from this figure that 

the decay width T(H — > ZZj) is several orders of mag- 
nitude smaller than those of H — > ZZ, H — > 77, and 
H — » jZ. For Higgs boson masses 190 GeV < m# < 
250 GeV, the ratios of the decay widths are T(H — » 
ZZj)/T(H -> 77) - T(H -> ZZj)/T{H -> jZ) < 10~ 7 
and T(H -> ZZ-y)/T(H -> ZZ) < 10" 11 . To iden- 
tify the origins of the smallness of these ratios, we note 
that, in addition to the cuts that we imposed on the de- 
cay products of H — > ZZ^j, which decrease the value of 
T{H — ► ZZj), there is also the suppression from three- 
body phase space, and from the higher order in the cou- 
pling constant a. These, however, do not completely ac- 
count for the suppression. There are other differences in 
the various decay amplitudes, which contribute to the 
small ratios. For instance, in the case of H — » 77, 
the decay amplitude receives contributions from charged 
fermion loops as well as a substantial contribution from 
W boson loops, whereas in the decay H — » there 
are no W boson loop contributions and the inclusion of 
anomalous triangle diagram, Fig.^c), further suppresses 
the amplitude. As a result of these differences, the sim- 
ple power counting method for estimating the size of the 
ratio of the decay widths T(H -> ZZ~f)/T(H -> 77) is 
rather unreliable. 

To investigate the dependence of the decay width 
T(H — > ZZj) on the helicities of the produced Z bosons 
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and the photon, we can use Bose symmetry and the CP 
invariance to obtain relations among the helicity ampli- 
tudes A\\'\^. Here, A and A' are the helicities of the Z 
bosons and A 7 is the helicity of the photon. As a con- 



sequence of these symmetries, 
satisfies the following relations 



the decay width T 



TaA'A^ = r_A_A'-A~ 



(4) 
(•5) 



In Fig.|21 we show the result of the calculation for the 
decay width T(H — » ZZj) as function of the Higgs boson 
mass rriH, for different helicities of the Z bosons and the 
photon. As it is clear from this figure, the decay widths 
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FIG. 3: (Color online) The decay widths T(H — > ZZ^f) as func- 
tion of ran for different helicities (AA'A 7 ) of the Z bosons and the 
photon are shown. The solid line is for the unpolarized case, the 

dashed line is for (0 + +), the dotted line is for (0 H ), and the 

dotdashed line is for (00+). The cuts imposed are the same as 
those for the total width T(H -» ZZ-y) of the Fig.|3| 
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FIG. 4: (Color online) The invariant mass distributions dT(H — » 
ZZ-y)/dm zz i as function of m zz i, the invariant mass of the final 
Z bosons, for Higgs masses of m H = 190, 195, 200, 210, 230, and 
250 GeV are shown. The cuts imposed are the same as those for 
the total width V(H -> ZZ-f ) of the Fig.0 
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FIG. 5: (Color online) The energy distributions dT{H —y 
ZZrf/dEj, as function of the photon energy Ej, for the Higgs 
masses of Fig.[I] are shown. The cuts imposed are the same as 
those for the total width T(H -> ZZ^) of the Fig.EJ 



To++ = To , which correspond to the case when one 

of the Z bosons is longitudinally polarized and the other 
Z boson has the same helicity as that of the photon, are 
the most dominant. This dominance is stronger for the 
higher Higgs boson masses. The decay widths = 

r — , r ++ _ = r — +, and r+_+ = r_+_ = r_ ++ = 

, that are not shown in the Fig.0 are negligible. 

This pattern of the polarization states may be viewed as 
a signature for the decay products of the process H — > 
ZZ^, since a longitudinal-transverse helicity combination 
cannot occur for the Z pair from H — > ZZ. 

In Fig.QJ we show the invariant mass distribution 
dT{H — > Z Z"f)/dmzz' as function of the Z boson pair 
invariant mass mzz 1 , and in Fig-El we show the energy 
distribution dT(H —y ZZ^)jdE 1 as function of the pho- 
ton energy E 1 . Notice that the study of the invariant 
mass and energy distributions for the Higgs boson masses 
of about 190 GeV is difficult. This is due to the imposed 
cuts that make the phase space quite small. 



3. SUMMARY AND CONCLUSIONS 

In the standard model, the three-body decay of the 
Higgs boson H — * ZZ7 is highly suppressed. However, 
this decay mode has some interesting features that sepa- 
rate it from decay modes such as H — > 77 and H — > jZ. 
One is the absence of W boson contributions in any of 
the loops in the Feynman diagrams for the H — > ZZj. 
Its amplitudes are dominated by top quark loops and 
therefore sensitive to top-Z couplings. Also, there is the 
presence of an anomalous vertex in the s-channel Z ex- 
change diagram of the Fig.^c), which might be studied 
were it not for the smallness of the decay width. 

Our explicit calculations show that the most dominant 
helicity combinations for H — > ZZj occur when one of 
the Z bosons is longitudinally polarized and the other Z 
boson and the photon have the same helicity. This is a 
result that was not apparent at the outset. With enough 
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statistics, this feature might be used to discriminate the 
Z pairs of H — » ZZ^ from those of H — > ZZ, since the 
helicities of the Z pair in the latter decay cannot be in 
the combination longitudinal-transverse. 

In summary, we find that the decay width for the pro- 
cess H — » ZZ7 is exceedingly small compared to those 
of H — > 77 and H — > 7Z, and that the suppression is 
greater than the simple phase space and coupling con- 
stant accounting might suggest. The experimental mea- 
surement of the decay width for H — > ZZ7 will be ex- 
tremely difficult, especially the identification of the helic- 
ities of the Z bosons and photon. Therefore, we expect 
not to detect a signal of this mode at the standard model 
level, and any detection is an indication of new physics 



beyond the standard model. 
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